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An approximate technique is presented for the evaluation of the 
mean and variance of the power sums with log-normal components. 
Exact expressions for the moments with two components are devel- 
oped and then used in a nested fashion to obtain the moments of the 
desired sum. The results indicate more accurate estimates of these 
quantities over a wider range of individual component variances 
than any previously reported procedure. Coupling our estimates with 
the Gaussian assumption for the power sum provides a characteri- 
zation of the cumulative distribution function which agrees remark- 
ably well with a Monte Carlo simulation in the 1 to 99 percent range 
of the variate. Simple polynomial expressions obtained for the mo- 
ments lead to an effective analytical tool for various system perform- 
ance studies. They allow quick and accurate calculation of quantities 
such as cochannel interference caused by shadowing in mobile te- 
lephony. 

I. INTRODUCTION 

The power sum with K independent components 

K 



P K = 10 log! 



X 10**" 



(1) 



is a random variable which appears in many areas of communications. 
With X k Gaussian, the quantity 

L k = 10** /10 (2) 

is called a log-normal variate. The characterization of the sum of L k is 
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of importance in multihop scatter systems, 1 log-normal shadowing 
environments, 2,3 target detection in clutter, 4,5 and the general problem 
of propagation through a turbulent medium.* Thus, the distribution 
and moments of Pk are quantities of considerable importance. Unfor- 
tunately, these quantities do not appear to be expressible in simple 
analytical formulae and, as a consequence, approximate procedures 
have been investigated for some time. Of particular interest is the 
Wilkinson approach which uses a normal approximation for the distri- 
bution of Pk. The problem of characterizing the distribution function 
then reduces to finding the first two moments of the power sum. 

The Wilkinson approach is consistent with an accumulated body of 
evidence indicating that, for the values of K that are of interest, the 
distribution of the sum of a finite number of log-normal random 
variables is well- approximated, at least to first-order, by another log- 
normal distribution. 1,9_12+ The central question, then, is how to estimate 
the mean and variance of the approximately Gaussian variate Pk, i.e., 
the power sum. Knowledge of the mean, m x , and standard deviation, 
a x , of the Gaussian random variable Xk leads to a complete specification 
of the log-normal variate L*, assuming, of course, no location param- 
eter for the log-normal. Yet, the procedure for then estimating the 
moments of the power sum variable 

Pk = 10 log, (L) = 10 log, (Li + L 2 + • • • + L K ) (3) 

based on the individual moments m x , a x , is by no means straightfor- 
ward. As discussed in the next section, the Wilkinson approach leads 
to useful results only for a limited range of small values of the dB 
spread a x . Unfortunately, this is not the a x range of most practical 
interest. 

The purpose of this paper is to outline a procedure that appears to 
be more accurate and to have a wider range of applicability than 
previously reported approximations. The new technique can be easily 
summarized: Analytical formulae are developed to compute the exact 
mean and variance for the power sum with two components, 

P2=101og(Li+L 2 ). (4) 

We then suppose P 2 is Gaussian, i.e., it is assumed that L\ + L 2 = 
10 P2/I ° is log-normal. We then consider 

P 3 = 10 log(L, + L 2 + L 3 ) = 10 log(10 P2/ '° + L 3 ) (5) 

and compute the mean and variance using the derived formulae. In 
this manner, we iterate until the required moments of Pk are obtained. 



* Many more applications, including the field of economics and the estimation of 
crude oil reserves, are noted in Refs. 6, 7, and 8. 
f We also confirm this observation here by means of Monte Carlo simulations. 
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We have compared the results of a Monte Carlo study to the 
estimates of moments obtained from the new method just described, 
and permuted versions of it. Our procedure leads to remarkably 
accurate estimates for the mean of Pk, over a wide range of component 
variances. The estimate of the standard deviation of Pk is accurate for 
a somewhat more limited range but, nevertheless, a much wider range 
than for previously reported approximations. 

The rest of this paper is organized as follows. In the next section, we 
establish additional notation and discuss, in more detail, previous work 
relating to power sums. The analytical development is presented in 
Section III, with the details reserved for the Appendix. The develop- 
ment of both accurate and simple expressions for the mean and 
variance of P2 is of considerable importance for system performance 
studies in areas such as shadowing in mobile radio and log-normal 
fading in microwave radio. In Section III, we also present such an 
expression in the form of simple polynomials. Simulation results, 
validating both the Gaussian approximation for Pk and the new 
estimation technique, are presented in Section IV. Concluding com- 
ments are given in Section V. 

II. PRELIMINARIES 

The power sum, as defined in (1), is measured in decibels. However, 
rather than use the definition for L* given in (2), it is more convenient 
to use the natural logarithm: 

L, = e Y *. (6) 

The relationship between the two associated normal variates is simply 

Y k = XX k , (7) 

where 

A = Vio logelO = 0.23026. (8) 

When the mean and variance of Xk are specified as, say, m x and a x , in 
order to use the representation (6), we have the obvious scaling 

m y = Xm x , Oy = X 2 a x . (9) 

Similarly, we define the quantities 

L = l0 p «/» = e z , (10) 

and the corresponding moments are related by 

m z = Xmp K , o\ = X 2 op k . (11) 

The rth moment of a log-normal variate L is given in terms of the 
moment-generating function of Y: 
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E(Li) = E(e rY >) = e"^,+V2(rV.>. (12) 

Note that the moments grow exponentially with the order r 2 . This 
rapid growth of moments may be one of the reasons for the limited 
range of applicability of the classical Wilkinson approximation (see 
Ref. 11); dealing with differences of large numbers can lead to numer- 
ical instability in trying to evaluate the approximate log-normal density 
of the sum. 

Marlow 13 has shown that under quite general conditions, as K — > 
00, P K is asymptotically normally distributed.* For fixed, finite K, and 
with the Xk independent, identically distributed normal random vari- 
ables (with mean m x and standard deviation a x ), Marlow also derived 
a small variance normal approximation. 14 As a x —* 0, a scaled version 
of Pk approaches the distribution function of a unit normal. 

In terms of log-normal variates, if we write 

P K = 10 logio(L ) = 10 logio(Li + L 2 + • • • + L K ) , (13) 

Marlow's second result says that, in the finite-component, small-vari- 
ance case, the sum of log-normal variates is also approximately log- 
normal. This result forms the analytical basis for the classical Wil- 
kinson approximation, which is the usual normal approximation for 
the distribution of the power sum. 
The Wilkinson approximation proceeds as follows: With 

L = e z , (14) 

now taken to be log-normal, Z is the associated Gaussian random 
variable with parameters m z and a z . To find these quantities, one uses 
(12) and equates the first two moments of both sides of (14). For 
example, with K = 2, and Y h and Y 2 identically distributed with mean 
m y and variance a 2 ,, 

E(L) = e m ' +m " i = 2e m > +1/2 ^ 

E(L 2 ) = e 2m ' +2ai = 2e 2m ' +2 °? + 2(e m > +i/2 ^) 2 . (15) 

Taking logarithms gives a set of linear equations for the two unknowns, 
m z and oi. 

Figure 3, which will be discussed in detail later, presents the results 
of a Monte Carlo simulation and the Wilkinson approximation, em- 
ploying eq. (15). As indicated, the Wilkinson approach tends to break 
down for o x = a y /X greater than 4 dB. This is consistent with the 
above-quoted result from Marlow. 



*For the infinite component case, Marlow also points out that there is asymptotic 
normality both on the power scale, as well as on the dB scale. That is, both Pk and L are 
asymptotically normal. 
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If one is interested in small a, but larger values of the random 
variable, i.e., the tails of the power sum distribution, Fenton 1 suggests 
a better approximation to the equivalent log-normal distribution by 
utilizing higher-order moments. In this case, eq. (15) is replaced with 
equality of, say, third and fourth moments. This procedure of using 
different values of r in different regions results in an approximating 
distribution which would be a series of connected straight lines when 
plotted on log-probability graph paper. 

At the other end of the a x component range, Farley 9 has derived a 
large variance approximation. With independent, identically distrib- 
uted variables, as a x — * °° 



Pr 



— (P K -m x )<a 

Or 



= m<x)] K , (i6) 



where <!>(«) is the unit normal distribution function. (Clearly, this is a 
decidedly non-Gaussian result.) Farley's simulations indicate that the 
large variance approximation gives better results than Wilkinson's 
approach for about a x > 10 dB and K fixed. 

In this study, we also take K finite, but focus on the midrange, 4 < 
a x < 12 dB, where neither Marlow's nor Farley's results are applicable. 
This midrange is of particular interest for the shadowing phenomena 
in mobile telephony, microwave radio fading, and airborne radar 
clutter. 

III. ANALYTICAL RESULTS 

We assume that 



K K 



L^L*^^^ (17) 



*=i 



for finite K is log-normal. The associated Gaussian variate is Z, and 
the goal is to determine the quantities 

m 2 = E(]n L) 

o\ = E[(ln L - m 2 ) 2 ] (18) 

in terms of the mean and standard deviation of the component variates 

Yk. 

Below, we give exact expressions for m z and o z when there are two 
components. The procedure for K greater than two combines the log- 
normal variates in a nested fashion two-by-two, using the exact for- 
mulas developed. To illustrate the procedure, suppose we have three 
log-normal variates: L = L\ + Lz + L3. The exact first and second 
moments of ln(Li + L2) are computed. We then take 

Z 2 = ln(L, + L 2 ) (19) 
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as a normal random variable and write* 

Z = ln(L) = M<? 2 + L 3 ) = ln(e Z2 + e Y »). (20) 

We then compute the moments of ln(L), again using the exact expres- 
sions developed for two components. This procedure has been tested 
in a variety of situations and groupings of variates. The results are 
very encouraging and will be discussed in the next section. Here, we 
outline the development of the analytical formulae, with the details in 
the Appendix. 

Naus 15 has derived the moment-generating function, and computed 
the first two moments of the power sum with two independent, 
identically distributed normal components. (Hamdan 16 generalized 
Naus' result for correlated variables with unequal variances.) Here, we 
extend Naus' general approach to the case of two independent vari- 
ables with unequal means and variances. 

For the two components, we have 

e 2 = L = Li + U = e v ' + e Y > (21) 

or 

Z = ln(e y ' + e y2 ). 

Define the Gaussian random variable 

w = Y 2 - Yi (22) 

m w — w = m.y 2 — m yi (23) 

<£ = <4 + oV (24) 

Taking the expectation of Z, we have 

E{Z) = £[ln(e Y| + e Y *)] 

= E{]n[e Y >(l + e y2 " y ')]} 

= E(Yi) + £[ln(l + e w )]. (25) 

The second term is 

£[ln(l + e w )] = I [ln(l + e w )Y(w)dw, (26) 

where f(w) denotes the normal density with the above-indicated mean 
and variance. The logarithmic term is now expanded in a power series. 
To ensure convergence of both the power series and the series resulting 
from the subsequent integration, the integral is broken up into the 
appropriate ranges. Thus, 



* Equivalently, L\ + L% is assumed to be log-normal. 
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I ln(l + e")f(w)dw = ln(l + e w )f(w)dw 



[ln(l + e~ w ) + w]f(w)dw. (27) 



The expansion 



ln(l + x) = X C jX \ Cj = ±—!— , (28) 

7=1 J 



valid for \x\ < 1, can now be used for each of the above integrals. The 
result (details are given in the Appendix) for the first moment is 

m z = E[]n(Lx + L 2 )] 

= my, = e- m " /2ai + m u — 






g*"..^"^"^ 



Ow 



m 



- ka 2 

W KUu 



Ow 

= m. V| + Gi(o M ,, m u ) , (29) 

where O(^) is the Gaussian probability distribution function as defined 
in the Appendix, eq. (39), and G\ is defined implicitly. It is not difficult 
to show that the series in (29) converges. Our computer experience 
indicates that about 40 terms are required for the 4th significant digit. 

As might be expected, the formula for the second moment is more 
complicated and is given by eq. (72) in the Appendix. The basic idea 
however, remains the same: When expansions of logarithmic functions 
are used, the expectation is broken up into integration ranges which 
provide convergent series. 

Although the formulae as developed involve infinite series and are 
complicated, their application is relatively simple. We can write [see 
eqs. (49) and (73) in the Appendix]: 

m 2 = m y , + Gi(a„, m„) 

oi = a?, - Gf(o„; m„) - 2p 2 G.i(a„., m„.) + G2(o w , m u ), (30) 
where the G,(o w , m„) are defined in the Appendix and 

p = -OyJOu,. (31) 

Observe that the G, functions, i = 1, 2, 3, depend only on the two 
parameters m w , a w . Consequently, we can evaluate Gi(a u ., m w ) as a 
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Table la — Coefficients of the approximating polynomials — Region I 
[-20 < m w < 0, < a w < 1 5] 







Functions 




Coefficients 


Gl 


G2 


G3 


AOO 


-0.1153239E 00 


0.4012876E-01 


-0.3958699E 01 


A01 


-0.5667912E 01 


-0.4483259E 01 


-0.5454983E 01 


A02 


0.1151279E 02 


0.7391760E 01 


0.1139280E 02 


A03 


-0.7162489E 01 


-0.3772190E 01 


-0.7116366E 01 


A04 


0.1312986E 01 


0.5262268E 00 


0.1315218E 01 


A10 


-0.1611385E 00 


-0.1579114E 01 


0.6839918E 01 


All 


0.20842 15E 02 


0.1637249E 02 


0.1962529E 02 


A12 


-0.4499768E 02 


-0.2994901E 02 


-0.4314091E 02 


A13 


0.27562 10E 02 


0.1337812E 02 


0.2647795E 02 


A14 


-0.5109783E 01 


-0.1747459E 01 


-0.4940592E 01 


A20 


0.1345124E 00 


0.2174588E 01 


-0.4717296E 01 


A21 


-0.2670183E 02 


-0.2114152E 02 


-0.2486802E 02 


A22 


0.5919191E 02 


0.3976749E 02 


0.5618572E 02 


A23 


-0.3695334E 02 


-0.1776890E 02 


-0.3501815E 02 


A24 


0.6912766E 01 


0.2195322E 01 


0.656O938E 01 


A30 


0.8057054E-01 


-0.7525302E 00 


0.1819360E 01 


A31 


0.1429709E 02 


0.1147157E 02 


0.1323538E 02 


A32 


-0.3225969E 02 


-0.2189501E 02 


-0.3050022E 02 


A33 


0.2055628E 02 


0.10O6767E 02 


0.1937757E 02 


A34 


-0.3888684E 01 


-0.1224983E 01 


-0.3668379E 01 


A40 


-0.3145306E-01 


0.8447987E-01 


-0.2817493E 00 


A41 


-0.2730047E 01 


-0.2222839E 01 


-0.2518272E 01 


A42 


0.6244253E 01 


0.4289557E 01 


0.5893337E 01 


A43 


-0.4048245E 01 


-0.2035757E 01 


-0.38O9492E 01 


A44 


0.7746786E 00 


0.2499568E 00 


0.7260144E 00 



function of these two parameters and store the results in three lookup 
tables. An alternative is suggested by Fig. 5, which gives the three G.'s 
as a function of a w with m w as a parameter. Since the curves are 
smooth, we can fit low-order polynomials to them and use the resulting 
expressions in (30) for the evaluation of the moments m z and a z . The 
advantage of utilizing these analytical expressions in system perform- 
ance studies is clear and, indeed, is what we have implemented. 

A least-squares fit was performed to determine the polynomial 
coefficients. After some experimentation, we chose the polynomial 

logioG.Ca,,,, m w ) = £ £ A Jk (i)oi /2 \m w \ k/2 , i = 1, 2, 3, (32) 

with J = K = 4. The coefficients Ajk(i) were obtained for two regions 
of the parameter space. The results for Region I (—20 < m w < 0, < 
a w ^ 15) and Region II (—40 < m w ^ —20, < a w < 15) are presented 
in Tables la and lb. Over these regions of interest, the maximum error 
produced by using the polynomial approximations was found to be 
about 1 percent as compared to the exact calculations based on eqs. 
(49) and (73). 

Before discussing the results of the simulation study, we digress to 
note how the analytical results can be extended. First, the development 
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Table lb — Coefficients of the approximating polynomials — Region 
[-40 <m w < -20, < a w < 1 5] 







Functions 




Coefficients 


Gl 


G2 


G3 


A00 


-0.3354792E 03 


0.8942737E 02 


-0.3596955E 03 


A01 


0.5281230E 03 


-0.1521976E 03 


0.5576596E 03 


A02 


-0.30854 14E 03 


0.9734654E 02 


-0.3242909E 03 


A03 


0.7895054E 02 


-0.2744724E 02 


0.8271772E 02 


ACM 


-0.7518427E 01 


0.2756644E 01 


-0.7863267E 01 


A10 


0.1354559E 04 


-0.3688427E 03 


0.1443269E 04 


All 


-0.2133838E 04 


0.6287175E 03 


-0.2250748E 04 


A12 


0.1 24549 IE 04 


-0.4077784E 03 


0. 13071 18E 04 


A13 


-0.3191052E 03 


0.1148104E03 


-0.3336339E 03 


A14 


0.3036738E 02 


-0.1161187E 02 


0.3167827E 02 


A20 


-0.1865732E04 


0.4966108E 03 


-0.1986007E 04 


A21 


0.294 1369E 04 


-0.8514775E 03 


0.3106047E 04 


A22 


-0.1718295E 04 


0.5574036E 03 


-0.1804377E 04 


A23 


0.4403 198E 03 


-0.1588548E 03 


0.46034 19E 03 


A24 


-0.4188256E 02 


0.1624075E 02 


-0.4366067E 02 


A30 


0.1059090E 04 


-0.2684747E 03 


0.1128759E 04 


A31 


-0.1671264E 04 


0.4638312E 03 


-0.1768423E 04 


A32 


0.9776851E 03 


-0.3061915E 03 


0.1028409E 04 


A33 


-0.2507 169E 03 


0.8842364E 02 


-0.2624284E 03 


A34 


0.2384698E 02 


-0.9152624E 01 


0.2487490E 02 


A40 


-0.2121206E 03 


0.5077676E 02 


-0.2264791E 03 


A41 


0.3350888E 03 


-0.8833020E 02 


0.3553536E 03 


A42 


-0.1963238E 03 


0.5876009E 02 


-0.2069294E 03 


A43 


0.5040234E 02 


-0.1717538E 02 


0.5284466E 02 


A44 


-0.4795980E 01 


0.1799387E 01 


-0.5008795E 01 



is valid when the random variables Y\ and Y2 are correlated as well as 
having different means and variances. With a simple modification in 
eq. (24) to account for the correlation [and another in (52)], 



07j, = O'y., + Oy t — PnOyfiy.,, 



(33) 



the results remain unchanged. Furthermore, in principle, the under- 
lying component variables do not have to be Gaussian. The proviso, of 
course, is that the integration with the non-Gaussian f{w) be tractable, 
and it must fall off fast enough so that the resulting series is convergent. 
(Clearly, the components are no longer log-normal.) 

Finally, we observe that the results can be extended in yet another 
direction. By generalizing the procedure outlined above, exact formulae 
for moments of more than two components can be obtained. However, 
the resulting expressions are very complicated and, given the high 
accuracy of our new method, the added complexity does not appear to 
be warranted. 



IV. SIMULATION RESULTS 

We use the definitions introduced earlier: 
U = 10** /10 



(34) 
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Pk = 10 logio(L) = 10 logio 



K 



k=l 



(35) 



In the following discussion, by m p and o>, we shall mean the first 
moment and standard deviation of the random variable Pk- The value 
of K should be clear from the context. We focused our Monte Carlo 
study on three issues: 

(i) How good is the assumption that L = L\ + L2 + • • • + Lk is 
approximately log-normal or, equivalently, how well is the cumulative 
distribution function (cdf) of Pk = logio(L) described by a Gaussian 
cdf ? Furthermore, how closely does the Gaussian cdf (based on the 
calculated m p and a p ) match the true cdf of the power sum? 

(ii) How accurate is the new analytical method in estimating the 
resultant mean m p and standard deviation a p and for what range of dB 
spread a x and component number K is the technique accurate? 

(Hi) Is the nested procedure we have described numerically robust, 
or are the estimates of m p and a p sensitive to the order in which we 
combine the log-normal components? 

The results of our Monte Carlo simulation confirmed the observation 
of a number of other investigators, that the cdf of P K is well approxi- 
mated by the Gaussian cdf, particularly in the range of practical 
interest.* Shown in Fig. 1 is the cdf of the sum of two log-normal 
variates with a x = 6, 10, and 14 dB. We note that the Monte Carlo 
simulation agrees quite closely with the assumed Gaussian cdf (based 
on the calculated m p and a p ) in the range of 0.1 to 99 percent. Outside 
this range, the simulation values start to deviate slightly. The cdf of 
the sum of a large number of log-normal variates, discussed below in 
Examples 1 to 3, together with a fourth case (the sum of six equal 
components with m x = dB and a x = 10 dB), are presented in Fig. 2. 
Here we observe excellent agreement between the calculated (as- 
sumed) and simulated cdf in the range 1 to 99 percent. 

Having confirmed the Gaussian approximation for Pk and the fact 
that the calculated cdf closely represents the true cdf, we turn our 
attention to the second issue, estimating the moments m p and a p from 
the moments of the individual log-normal components. In Figs. 3a to 
3c, the number of identically distributed components, K, is held fixed 
and the component variance (dB spread) a x is varied. In all cases, the 
component mean m x is 0. Figure 3a illustrates the results for two 
components and serves to verify the Monte Carlo program. Our 
computed results coincide (as they should) with the simulation. In 
contrast, the Wilkinson approximation begins to give inaccurate results 



* The simulations typically had 10,000 sample points for each value quoted or 
indicated on the figures. 
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99.8 



99.99 



Fig. 1 — Cumulative distribution function of sum of two log-normal variates with zero 
means and variances equal to 6, 10, and 14 dB. 
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1. THREE VARIATES: m = 0, ff = 6, 7,9,5. 

2. SIX EQUAL COMPONENTS: m - 0, CT- 10. 

3. NINE COMPONENTS IN THREE GROUPS: GROUP 1,/n - -10. O* 6; GROUP 2, m = -18. 

0-= 10; GROUP 3. m - -38, <7=" 12. 

4. EIGHTEEN COMPONENTS IN THREE GROUPS: GROUP 1, m = 10, CT- 10; GROUP 2, m - -2, 

<T= 10; GROUP 3, m ■ -8, cr- 10. 

Fig. 2 — Cumulative distribution function of the sum of a large number of log-normal 
variates corresponding to Examples 1 to 3 in Section IV plus an additional case. 
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at about a x = 4 dB and becomes increasingly worse with increasing 
a x . 

Figure 3b illustrates six components, and Fig. 3c illustrates ten. The 
dB spread ranges from 2 to 14 dB. Our technique gives remarkably 
accurate estimates of m p , while the estimate of a p is quite close for 
o x < 8 dB and then tends to underestimate the true standard deviation. 
By way of contrast, for these ranges and number of components, the 
Wilkinson approach gives gross inaccuracies and cannot be used. 

In Figs. 4a and 4b, we vary the number of components, set m x = 0, 
and hold the dB spread fixed. Again, our procedure gives very accurate 
estimates of the mean and, as the number of components increases, 
the variance is underestimated only slightly for a x = 6 dB (Fig. 4a) and 
somewhat more so when a x = 10 dB (Fig. 4b). For these cases, the 
error in the estimate of the variance as a function of the number of 
components is as shown in Table II. (There are negligible errors in 
estimates of the mean.) As indicated in the figures, the errors in the 
Wilkinson approach are substantial, ranging up to 74 percent when 
a x = 10 dB. 

It should be observed that the situation with equal mean components 
and large dB spread is probably a worst-case situation. Indeed, from 
the above table we see for K = 16, a x = 10 dB, that the estimate of a p 
is off by 12.8 percent. In contrast, in Example 3 below, we deal with 
K — 18 components, and the same dB spread. There, the error in the 
estimate is less — about 6 percent. The only substantial difference in 
the two situations is the unequal mean values. These comments 
notwithstanding, further study is warranted for this boundary of 
parameter values, i.e., where the number of equal mean components 
is high (K > 16) and the dB spread is identical and large (a x > 10 dB). 

Our technique was also tested on a number of examples with unequal 
components. We report on three that were chosen to illustrate the 
general applicability of the new method. 
Example 1 — Three components, equal means, different variances: 

m Xi = 0dB, i-1, 2, 3; a Xj = 6, a X2 = 7, a l3 = 9.5dB. 

The Monte Carlo simulation gave m p = 8.08, a p = 5.356. Our procedure 
provided the estimates m p = 8.05, a p = 5.273, for a 0.03-dB error in the 
mean and 1.5 percent error in the variance. 
Example 2 — Nine components, taken in homogeneous groups of three: 



m Xi = -38 dB 
m,, - -18 dB 
m* ; = -10 dB 



a Xi = 12 dB, i = 1, 2, 3 
a Xi = 10 dB, i = 4, 5, 6 

o I . = 6dB 1 i = 7, 8, 9. 
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Fig. 4— Mean and standard deviation of power sum with varying number of identically 
distributed components, (a) Component mean, m x = and standard deviation, a x = 6 
dB. (b) Component standard deviation a x = 10 dB. 



Monte Carlo results gave m p = -0.61, a p = 3.90, while our technique 
yielded m p = -0.6, a p = 3.79. This gives an error of 0.01 dB in m p and 
2.8 percent in o>. 

We now address the third issue raised at the beginning of this 
section. To show that the way in which we combine the variables has 
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Table II- 


-Error in estimate of variance 








No. of Components (K) 




Percent Error 


2 


4 8 16 


32 


o, = 10 dB 
a, = 6 dB 


0.13 
0.43 


1.4 0.4 12.8 
2.0 3.3 8.2 


21.4 
11.5 



little effect, we permuted the order in combining the nine components. 
The resultant mean estimate varied between -0.59 and -0.64 dB, 
while the standard deviation varied between 3.66 and 3.89 dB. This 
example — and a number of others studied — illustrates the desired 
numerical robustness of the procedure. 

Example 3 — Eighteen components: This example would correspond to 
two tiers of interferers in a cellular mobile radio scheme. 

Components 1 to 6: m x = 10 dB, a x = 10 dB. 
Components 7 to 12: m x = -2 dB, o x = 10 dB. 
Components 13 to 18: m, = -8 dB, o x = 10 dB. 

The Monte Carlo results are: m p = 27.07, o p = 4.54, and our technique 
gave m p = 27.04, a p = 4.26 . The error is 0.03 dB in the mean estimate 
and 6.2 percent for o p . 

V. CONCLUSIONS 

In this study, we have verified once again, as others have, that the 
sum of a moderate number of log-normal random variables is well 
approximated by another log-normal variate, especially in the cdf 
range of 1 to 99 percent. Perhaps more important, we have been able 
to relate the mean and variance of the resultant (Gaussian) power sum 
to the first two moments of the individual underlying Gaussian com- 
ponents. Our method is highly accurate in the range of parameter 
values of most practical interest. 

Based on our preliminary simulations, we may conclude that the 
analytical method presented in this paper is extremely accurate in 
evaluating the mean of the power sum for the complete range of 
parameter values investigated. This range was a dB spread of 2 < 
a x < 14 dB and up to K — 30 components. It is also accurate in 
estimating the standard deviation with up to eight equal log- normal 
variates when a x is less than 10 dB. Outside this range, the approxi- 
mation is slightly less accurate. 

APPENDIX 

Derivation of First and Second Moments 

Our purpose is to provide enough detail so that an interested reader 
can rederive the formulae for the first and second moments, as given 
below by eqs. (48) and (73). 
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The expansions we use are 

ln(l + *) = £ C***, |*|<1 and C k = , (36) 

A-l * 

ln 2 (l + x) = J 6*x**\ 
|*| <1 and b k = \' IT 1 . (37) 

R + I y— 1 

The density function of a normal random variable with arbitrary mean 
and variance is defined by 

*' (*; m, a) = -^= e~ (x ~ m )2/2a \ (38) 

V27ra 

The distribution function is a simple parameter function and is given 
by 



(39) 



The function <&(x) is related to the error function by 

erf(*)=4= F e ~* dt i40) 

vw Jo 

O(x) = x h + V4 erf(x/>/2). (41) 

Three integrals which occur often in the development are 



/•+00 



i 



E{e rx ) m e rx &(x; m, a)dx = e mr+r ° /2 (42) 



° e rx <t>'(x;m,a)dx=[e mr+r ^ /2 ]^ m ^ 



I 



H e- rx V(x; m, a)dx = [e^^W^-^ . (43) 



We repeat the definitions introduced in Section HI: 
e 2 = L = Li + L 2 m e Yl + e Y * 
Z = ln(e Y ' + e Yi ), 
with the Yi Gaussian random variables. Let 
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w=Y 2 -Y l 

m w = w = E(w) = m.y 2 — m yi 

al = E(w - w) 2 = o& + o&. (44) 

For convenience, we shall let f(w) denote the particular Gaussian 
density function 

f{w) = $>'{w; m w , a w ); (45) 

Z is rewritten as 

Z = ln(e y ' + e Y *) = ln[e y '(l + el]- (46) 

Taking the required expectations 

E{Z) =m z = m yi + E ln(l + e w ) 

= m yi + I ln(l + e w )f(w)dw. (47) 

As discussed earlier, the integral is broken up into the intervals 
(—oo, 0) and (0, oo). Upon application of the series (36) and the 
averaging (42) and (43), we obtain the result 

m z = rriy, + m w ^\ — ■== e- m - /2a - 

\°u>J V27T 



+ 2 C*e*^ 



— m w — kof 



m w — kol 



+ e - km «<& 

All the terms, except m, lt are lumped together and denoted by 
G\(o w , m w ): 

m z = m yt + Gi(a w , m w ). (49) 

The function Gi is shown graphically in Fig. 5. 
The second moment is decomposed as 

E(Z 2 )=E{[Y 1 + ln(l + e")] 2 } 

= o 2 + m 2 , + £[2Y, ln(l + e w )] + £[ln 2 (l + e w )]. (50) 

First, consider the cross term 

A = E[2Yi ln(l + e w )], (51) 

where Y\ and w = Y2 — Y\ are jointly Gaussian with correlation 
coefficient 
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Fig. 5 — The plots of G,'s as a function of a w with m w as a parameter. 
(Yi - m y ){w - m w ) 



p = E 



Jy^w 



= -Oyjy/ol + 



y, v u_y, -r vy 2 



(52) 



We utilize the conditional expectation property of jointly Gaussian 
random variables (see Ref. 17, Section 7-5): 



E[(Yi - m y )\{w — m w )] = — — (w — m w ) 



(53) 



or 



E[(Yi\w)] = m yt - p 2 (w - m w ). 
Then, (51) becomes 
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A = E[2Y 1 ]n(l + e w )] 

= 2E w []n(l + eT)E Yl \w(Yi)] 

= 2E w {]n(l + e w )[m yi - p 2 (w - m w )]} . (54) 

For convenience, we denote the first and second terms by A\ and A 2 , 
respectively: 

A=Ax-A 2 . (55) 

Since 

Ai = 2m >1 £[ln(l + e w )] 

= 2m yi (m z - m yi ) (56) 

from (47), we need only consider the second term A2, 

A 2 = 2p 2 E[(w - m w ) ln(l + e w )] . (57) 

Integrate A 2 by parts: 



= 2p^ 



(w - m w ) ln(l + e w )f(w)dw 



= 2p 2 o 2 A T ^f(w)dw. (58) 

For negative w, we use the expansion 

= £ (-l) k e tk+1)w (59) 



1 + e 

to obtain 



J 1 + e w 



At. = 2p 2 al \ -f(w)dw 



= 2p 2 al £ (-i)V"- ( * +1)+( * +1,1 °- /2 



<E> 



A-0 
— JTlw — 



AC* + 1) 1 

a w I 



For the other part of (58), we have 



e 
Ai" = 2p 2 o 2 w I ——f(w)dw 

)o l + e 



I'M. 



(60) 



= 2p 2 al I , ; _ w f(w)dw 
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= 2p 2 o 2 | 2 (-D k e- kw (f(w)dw 

k=0 



= 2p 2 al J (-1)V*^ + *^ 2 






(61) 



Collecting the results for the cross term, 

A - Ax - A 2 = Ax - (A' 2 + A?) 
= Ai - 2p 2 G 3 {o w , m w ) 
= 2m yi G\(a w , m w ) — 2p 2 Gz(o W) m w ). 



(62) 



The remaining term in E(Z 2 ) to evaluate is the last expression in (50). 
We denote it by G 2 (o w , m w ). 

G 2 (a w , m w ) = E[ln 2 (l + e w )] 



ln 2 (l + e w )f(w)dw + ln 2 (l + e w )f{w)dw 
Jo 



= Bi + B 2 . 



(63) 



For the integration over the negative real line, we can use the expansion 
given by (37): 



B,= £ b k \ e {k+l)w f(w)dw 



= Y 6 A el<* +1 > m ">+(*+i> 2 °«/2] 
*=1 



<D 



— m w — (^(k + 1) 



O u - 



(64) 



The expression B 2 is rewritten as: 

li, | {]n[e w (l + e- w )]} 2 f(w)dw 



■f 



[ln 2 (l + e~ w ) + 2w ln(l + e~ w ) + w 2 ]f(w)dw 



= B*, + B\ 4- 53. 
The expressions for B3 and B5 are relatively straightforward: 



(65) 
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B 3 = w 2 f(w)dw 
'o 



= TWl 



1-0 



— m L 



+ V2/77 m w a w e- mi/2oi 



+ 2 



l-O 



-mi 



—TJlwO, 



Ou> / V^ 






B 5 = ln 2 (l + e- w )f(w)dw 



ln 2 (l + e w )f(-w)dw 



£ b k e ik+l)w f(-w)dw 



■I 



m w — ai(A + 1) 







The final term to consider is 

B 4 = 2u; ln(l + e- u ')f(w)dw 



k=\ 



= -2 £ C* we kw f(-w)dw. 



(66) 



(67) 



(68) 



Complete the square for the expression e kw f(—w) and integrate to 
obtain 



B 4 = -2 X c*e ( - m «* + * Jo " /2 » 



m*0 



-m* 



>/2tt 



g-mt/2o' 



, (69) 



where 

rrik = —m w + kcrl,. (70) 

Collecting the terms that make up G 2 (o w , m w ) in (63), we have 

G 2 (a w , m u ) =B i +B 2 = B l + (B 3 + B 4 + B 5 ) (71) 

as given by eqs. (64), (66), (69), and (67). 

To summarize, the second moment is given by the expression 
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E(Z 2 ) = 0^, + ml x + 2 m yj Gi(a w , m w ) 

- 2p 2 G 3 (o w , m w ) + G 2 (a w , m w ) , (72) 

and, finally, the variance is 
oi = E(Z 2 ) - ml 

= of, - G 2 (a w , m w ) — 2p 2 G 3 (a w , m w ) + G 2 (a w , m w ). (73) 
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